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Abstract 

We show the existence of quantum states of the Heisenberg XY chain which 
closely follow the motion of the corresponding semi-classical ones, and whose 
evolution resemble the propagation of a shock wave in a fluid. These states 
are exact solutions of the Schrodinger equation of the XY model and their 
classical counterpart are simply domain walls or soliton-like solutions. 
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I. INTRODUCTION 



The relationship between classical and quantum dynamics, although largely investigated 
since the beginning of quantum mechanics, is still a non exhausted problem of ever continuing 
interest. This has lead, for example, to the introduction of new concepts such as coherent 
states and their generalizations, which has been found to be successful in describing the 
quasiclassical behaviors of a wide class of physical problems. The possibility to construct 
quantum states which closely resemble the evolution of classical ones is appealing also in 
connection with the new developing area of quantum computing [JTJ] . In this context the basic 
component of a quantum computer is a q-bit i.e. a special quantum state which has to be 
manipulated by the logical quantum gates (unitary operators) and which has to be preserved 
'intact' for an extended period of time. In any practical task, however, the interaction of 
the system with the environment (heat reservoir) is unavoidable and decoherence usually 
develops. On the other hand, one may expect that, due to their intrinsic collective character, 
quantum states close to macroscopic (classical) ones are more robust against decoherence 
and therefore more suitable to store information in a quantum computer. The problem of 
finding such states, however, is far from being trivial even for isolated systems. Indeed, the 
complexity of a quantum system related to the huge dimensionality of the Hilbert space, 
represents a major obstacle in solving the time-dependent Schrodinger equation (one usually 
resorts to approximations which are not always well controlled). 

The aim of the present paper is to show that for the XY Heisenberg spin chain |2| it is 
possible to find quantum states which closely follow the motion of the corresponding classical 
ones and whose evolution resemble the propagation of a shock wave in a fluid. The classical 
counterpart of these states are moving domain walls or soliton-like solutions. Stable shock- 
wave solutions of different shapes were also encountered in the semiclassical dynamics of the 
Heisenberg chain || as well as in other similar models |4j]. Here we show, on the example 
of the solvable XY model, that a spin system may indeed have exact quantum states which 
recall classical moving shock-wave solutions. In particular, a detailed comparison between 
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the quantum and the classical evolution for several initial conditions is presented. The 
regions where the exact quantum states are well approximated by the classical ones are also 
investigated. In these regions the classical approach can be used to track the dynamics 
of those quantum observables which are technically unaccessible by a quantum-mechanical 
analysis. The paper is organized as follows. First we derive semiclassical and quantum 
equations of motion, and in particular solve them in the quantum case. This enables us to 
compare exact quantum and semiclassical time evolutions for a chosen set of initial conditions 
leading to shock-wave formation. The properties of the latter are then discussed, and in the 
conclusion the main results of the paper are summarized. 



II. SEMICLASSICAL ANALYSIS OF THE XY MODEL 

The quantum Heisenberg XY model [] is introduced as 

^ = ^E(« + i + « + i) (i) 

n 

where s l = \o % ■, % = 1,2,3 are spin 1/2 operators, a 1 are Pauli matrices, and J is 
the exchange constant related to the interaction between spins. As well known, the two- 
dimensional character of the interaction and its restriction to neighbouring spins leads to the 
exact solvability of this model || . In the following we shall derive semiclassical equations of 
motion of the XY model, by averaging the Heisenberg equation of motion ih^ = [/, H] for 
a generic observable / over (uncorrelated) coherent states. For s z n (site magnetization) and 
= s x ± is y (rising and lowering operators), we have: 
ds z J 

^~fu~ = ~r> ( S " S n+1 ~ S n S n+l ~ S n-l S n + S n-l S n) ) (2) 



= ±J< + , (3) 



x An arbitrary constant magnetic field BJ2 n s^ can be included; here we consider B = for 
simplicity. 



(in the following to simplify notation we fix % = 1 and J = 1). Introducing spin ~ coherent 
states as 

|A(t) >=n®M*) >, W{t) >= i ; = 1 T), (4) 
with | j) n denoting the spin up state at site n, we compute 
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where (s) denotes the average over the state |A). It is worth noting that these equations 
define an inverse stereographic mapping from the variables /i n , /z* , to the vectors S n (classical 
spins) on the unit sphere 

qx i^n ) {^n ) ny (^n)~j~ i^n ) qz / z \ (/?\ 

D n — 2 ' n ~~ 2 ' ~~ 

By averaging Eqs. @, © with respect to |A) and by eliminating redundancy (note that 
the quantities in Eq. (||) are not all independent), we finally obtain the equation of motion 
for fi n (t) as 



dt 2\ l + |/i n+1 | 2 l + |/i n _i| 2 /' 
Note that this equation (and its complex conjugated) can be put in hamiltonian form 



(7) 



fln = {Vn,H} (8) 

with respect to the Poisson bracket 

ft \ -S^fi i i \^\2 ( 9 f 9 9 df dg\ 



and the hamiltonian 



C ^'^ J 2 V(l + l/in| 2 )(l+|/in + l| 2 ) 1 ' 

These results were also derived in Ref. using the stationary phase approximation and 
the path integral formulation of quantum mechanics. It is remarkable that the classical 
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hamiltonian in Eq. ([H]), once expressed in terms of the vectors S n through Eqs. (|5]|J), takes 
exactly the same form as in Eq. (|]) with quantum operators replaced by classical spins. This 
is precisely what one would expect from the classical limit in coherent state representation. 
On the other hand, since coherent states do not preserve the hidden symmetry algebra 
of the quantum XY chain (they are not eigenstates of the monodromy operator) ||, the 
integr ability structure of the system is lost in the semiclassical limit. 

III. EXACT QUANTUM ANALYSIS OF THE XY MODEL 

In this section we shall derive exact analytical expressions for the time evolution of the 
expectation values (/(£)) = {^o l elW< /(0)e _i7 ^ f l^o) of observables f(t) of the XY model. To 
this end we take as initial state the totally uncorrelated quantum state 



hh>>= n 



/ e i<f„/2 cos (a n ) * 
ye-^sinK) J 



(11) 



By the Jordan- Wigner |7[ transformation from spin to Fermi operators 

s n. = II (l — 2C^C m ^ C n , s n = (s„)*, s z n = — — C^Cn, (12) 
the XY Hamiltonian ([!]) reduces to 

K = lY,( C n + lCn + C+C n+1 ). (13) 
z n 

This Hamiltonian is readily diagonalized as 

H = ^cos(k)4c k , (14) 

k 

by the Fourier transform C n = -^Y,k e ~ lkUc k- The Heisenberg equations of motion 

i ~dt = C0S (^ Cfc ' ( 15 ) 
are then solved as Cfc(t) = e - ttcos ( k ) Cfc(O) . Returning to the operators C n , we get 

C n {t) = J=^ e -*™^ CA; (0) =Y,i m - n Jrn-n(-t)C m (0), (16) 
V k m 



with J n (t) denoting the Bessel function of order n. From Eq. ([12]) and from its inverse 

Cn = ( II 4) 4, C+ = S-(U sA , (17) 

\m<n / \m<n / 

one can finally derive exact expressions for the time evolution of the observables of the 
original spin chain. Thus, for example, the time evolution of the averaged spin-z projection 
is obtained as 

s z P (t) = 2( s ;) = mi - 2c;(t)c P m ) = 

oo oo / m— 1 \ / n— 1 \ 

1-2 e e i n - m Jn- P (-t)j m - P (-t) this-! n n 4 

n=-oo m=-oo yi=-oo / y'=— oo / 

S p( q ) ~ E ^n(*) Jm(t)sin(2a fH ^)ain(2a g+n )oos f -(m -n) + $ m+P - $n+ P J JJ ( 1? 

(in a similar manner one proceeds for the other observables). In the next section we shall 
compare this expression for the average site magnetization, with direct numerical integra- 
tions of Eq. (0). 

IV. COMPARISON BETWEEN CLASSICAL AND QUANTUM EVOLUTIONS 

To compare the classical and quantum time evolution from an identical initial state, we 
shall first derive the classical initial conditions which correspond to the quantum state Eq. 



H|). To this end we remark that the expectation values of the operators s„ , s*, with respect 



to \ipo) are easily calculated as 

( s n) = ~^—sin 2a n , (s z n ) = -cos 2a n . (19) 

The initial conditions in the classical system are fixed by requiring that the above expressions 
for (s„), (s n ), coincide with the corresponding classical ones in Eq. @, this leading to 

/z n (0) = tgia^e-^ (20) 

as initial conditions for Eq.|7|. Since the classical states stay uncorrelated for all times while 
quantum ones develop correlations in the course of time, one expects the corresponding 



classical and quantum time evolution to be different for generic initial conditions. For special 
choices of the initial conditions however, these dynamics may qualitatively agree for times 
long enough to observe phenomena of the classical system, such as shock waves formation, 
also in the quantum one. Shock waves were observed in a semiclassical description of the 
Heisenberg chain using as initial conditions non constant (bell shaped) profiles for (s*) with 
a constant phase $ n = const along the chain ||. Here it is convenient to take as initial 
conditions constant values for the site magnetizations but not for the phase i.e. we use 

N N 

S*(0) = cos(2a), $„ = $(t0, ~-<n<- (21) 

with a an arbitrary constant and N the length of the chain. Note that, except for the case 
of linearly increasing phases, $ n = (fin (0 =constant) for which <S* will stay constant at 
all times (this follows directly from Eqs. (|7|jl8|)), the inhomogenuity of the phase at t = 
will in general induce a spatial dependence of S z in the course of time. In the following we 
choose the phase of the classical and quantum initial state so as to generate two types of 
excitations: localized states and expanding shock wave solutions. 

A) Localized states. 

We have generated these excitations with the following initial conditions 

a) = (fin, n < 0; = (fin + A, n > 0; ( dislocation), 

b) $ n = (fin, $ = A ^ 0; (local inhomogenuity). 

In Fig. p] we depict the time evolution of the quantum averaged site magnetization as 
computed from Eq. fll8|) for the state in Eq. (p[) with phase given by a). In Fig|| a 
direct numerical integration of Eq. (Jfty with a fourth order Runge-Kutta method for the 
corresponding initial condition in Eq. ( p0[) and $ n given by a), is reported. We see that in 
both cases localized soliton-like excitations are formed together with a radiative field which 
is much stronger in the case of the classical evolution. The initial conditions of type b) also 
lead to the formation of soliton-like solutions and background radiation, but in this case the 
quantum solitons have much shorter life-times and therefore, will not be discussed here. 

B) Expanding shock-wave solution. 
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To form a shock-wave in the magnetization profile we have considered initial conditions of 
the form $ n = 4>\n\, on an infinite open chain (note that the phase has a cusp at n — 0). Such 
a configuration can in principle be prepared by applying a periodic magnetic field cos (cut) 
to the left half of the chain and a corresponding 7r-shifted field (i.e. —cosiiot)) to the right 
half. In Figs. ^ |3] we depict the time evolution of this initial condition for, respectively, the 
quantum and the classical cases. The formation of an expanding symmetrical shock front is 
clearly seen in both cases and the corresponding time propagations are in a good agreement 
for long times. We also remark that the expanding inhomogeneous region with higher S* is 
caused by the fact that there are two opposite fluxes of magnetization j > (on the left) and 
j < (on the right), colliding at the center. At the initial time t = 0, they are separated by 
the single site n = 0. As time grows, the opposite magnetization currents become separated 
by the zero-flux region j « 0, expanding alongside with perturbed S* region as one can see 
from Fig.[5| The symmetry of the S* is just due to the symmetry of the incoming currents. 

The relation between the rate of change of magnetization in time and its current is given 
by a discrete continuity equation 
dSl 



n 



On) - O'n+l), jn = Z (4-l s n ~ s n-l s n) > (22) 



dt 

which is obtained from Eq. (0) by averaging with respect to the initial state. Denoting by 
—A(t),A(t) the boundaries of the perturbed region, and summing Eq. (|2^) over the points 
inside this region, we obtain 

= ~ On+l)) = (j-A(t)) - O'A(t))- (23) 

01 -A(t) 

For the case $ n = 4>\n\ we find 

~7 E <S£=(l- (^(0)) 2 ) sin{4>), (24) 

-A{t) 

from which we see that the excess magnetization in the perturbed region grows linearly 
in time 

A(t) 

E (Sn(t) - S*(0)) = at, a = (1 - (^(O) 2 )) S m(0). (25) 

-A(t) 
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Using a similar argument for the local energy 

E n = + ( 26 ) 

we obtain that the excess energy inside the perturbed region also increases linearly in time 

A(t) l 

£ (E n (t) - E z n (0)) =pt, (3 = — (l - (^(0)) 2 ) S 2 (O) S m(20). (27) 

-A(t) 

This is also seen from Fig.^| where the energy profiles of a shock solution at different instants 
of time are reported. It is worth to note that the magnetization production is maximal (see 
Eq. (j25D with = vr/2), when the energy production is zero (see Eq. ( p7|) ). So, quite 
surprisingly, the maximal magnetization current is always accompanied by a zero energy 
current. We observed this feature also in quantum evolution of the XXZ spin hamiltonian, 
for which the time-dependent Schrodinger problem can be exactly solved for special initial 
conditions ||. Another interesting feature of the above solution is its scaling property 

S z n (t) -> S z {n/t), n,t>l (28) 

for the site magnetization and energy, as well as for the flux and the energy j n (t), E n (t), 
which can be proved along the similar lines as in ||. Precise form of the scaling function 
is complicated and is omitted here. Finally, we note that Eq. (|28|) agrees with Eqs. (p5|). 
(fZTj), and that the above scaling properties are shared ( numerically, at least) by the classical 
solution. 

V. SUMMARY 

In this paper we have compared the exact time evolution of the quantum Heisenberg XY 
spin 1/2 model with its classical counterpart. In particular, we have reported on localized 
states and on expanding shock-wave solutions whose classical and quantum evolutions agree 
for long times. The problem of finding suitable conditions which allow to discriminate the 
quantum states which will evolve 'quasiclassically' from the ones which will not, remains 
however unsolved. Thus, for example, we found that initial conditions leading to decrease 
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of absolute value of local magnetization (e.g. initial conditions of Fig. f| with the sign of 
<p inverted), develop fast instabilities in the classical case which are not observed in the 
quantum evolution. This could possibly be explained by a difference of the 'dynamical 



temperature' for the classical spin system in the two cases 110]. It will be interesting to 



extend this analysis to other initial conditions such as the one used in Ref. [Ill], as well as 
to other quantum systems such as the XXZ model. 
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VII. FIGURE CAPTIONS 



Fig.l. Snapshots of site spin-z component s z n expectation values for localized excitation 
at a different moments of time, plotted from Eq. The initial conditions are: S z (0) = 

0.8, $ n = <pn,n < 0, $ n = <pn + 20, n > 0; 4> = tt/2. Graphs a,b,c,d correspond to 
time=10, 30, 50, 100 respectively. 

Fig. 2. the same as in Fig.l but for the semiclassical evolution obtained from Eq. (^). 
Graphs b, c, d correspond to time=30, 50, 100 respectively. 

Fig. 3. Average magnetization profile (s^) for shock solution (exact) at different times, 
computed from Eq. (|18"1). Consecutive expanding curves correspond to t = 10,30,50, 100. 
The initial conditions are: S z (0) = 0.8, $ n = <f)\n\, <fi = 7r/4. 

Fig. 4. the same as in Fig. 3 but for the semiclassical evolution. 

Fig. 5. Flux profile (exact) at time t = 0, 30, 50, 100. Initial conditions are the same as 
in Fig. 3. 

Fig. 6. Energy profile for a shock solution ( t = 0,30,50, 100). Initial conditions are the 
same as in Fig. 3. 
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